Relativistic Doppler effect: universal spectra and zeptosecond pulses. 
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We report on a numerical observation of the train of zeptosecond pulses produced by reflection 
of a relativistically intense femtosecond laser pulse from the oscillating boundary of an overdense 
plasma because of the Doppler effect. These pulses promise to become a unique experimental and 
technological tool since their length is of the order of the Bohr radius and the intensity is extremely 
high (X 10^^ W/cm^. We present the physical mechanism, analytical theory, and direct particle-in- 
cell simulations. We show that the harmonic spectrum is universal: the intensity of nth harmonic 
scales as 1/n'' for n < 47^, where 7 is the largest 7-factor of the electron fluid boundary, p = 3 and 
p = 5/2 for the broadband and quasimonochromatic laser pulses respectively. 
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Recent impressive progress in the physics of attosecond 
X-ray pulses Jj triggers a fascinating question whether 
a range of even shorter pulses is achievable with the 
contemporary or just coming experimental technology. 
In the present work, we show that when an ultra- 
relativistically intense laser pulse, /Ag > 10^°W/cm~^ 
/im^, is incident onto a plasma- vacuum boundary, the 
reflected radiation naturally contains subattosecond, or 
zeptosecond (zepto= 10~^^) pulses. 

The idea of harmonic generation via reflection from 
a laser-driven oscillating plasma boundary was first pro- 
posed a few years ago and has been studied in several the- 
oretical articles Numerous experimental works 
[EIEQ and simulations were devoted to the laser in- 
teractions with overdense plasmas and the high-harmonic 
generation. Plaja et al 3] were the first to realize that the 
simple concept of an oscillating plasma mirror gives an 
opportunity to produce extremely short pulses. In a re- 
cent work, Naumova et a/_4j have made a revolutionary 
step in this direction since they have proven that a single 
attosecond pulse can be isolated when a laser pulse fo- 
cused down to the A"^ volume is reflected from a plasma 
surface. 

The emphasis of our work is different. We concentrate 
on the very basic physics of the high harmonic genera- 
tion and discover the spectrum universality as well as its 
scaling. These features have been overseen in the earlier 
publications. The central result of our present work is 
the universal power-law harmonics spectrum decaying as 
l/riP. The exponent p — 5/2 for a quasi-monocromatic 
and p = Z for a broad band incident laser pulse. This 
power-law spectrum runs up to a critical harmonic num- 
ber Uc oc 47^^^, where an exponential cut-off sets in. 
Here 7inax is the largest relativistic factor of the plasma 
boundary. It is the slow power-law decay, that allows 
for the subattosecond pulses. We compare the analytic 
theory with direct particle-in-ccU (PIC) simulations. 

Let us consider a monochromatic laser wave with 
the dimensionless vector potential eAi{t, x)/mc^ = 
Re [ao exp(— iwt — icjx/c)]. This wave is reflected by a 
sharp plasma surface positioned at X{t') at the time t' . 



First of all, we must decide what boundary conditions 
should be used. We are interested in reflection from the 
overdense plasma N^, 3> Nc with F = aq (Nc/Ne) <C 1, 
where Ng is the plasma electron density, and Nc — 
(jj^me/47re^ is the critical density. 

The "ideal mirror" boundary condition implies zero 
tangential components of the vector potential at the mir- 
ror surface. As a consequence, when the ideal mirror 
moves with 7 3> 1 toward a laser pulse with the elec- 
tric field El and duration t, then the reflected pulse ac- 
quires the electric field Sj-efl °^ I'^Ei and the duration 
Tj.gg cx •r/7^. Consequently, the energy of the reflected 
pulse must be 7^ times higher than that of the incident 
one. However, as the plasma surface is driven by the 
same laser pulse, this scaling is energetically prohibited, 
and the plasma cannot serve as an "ideal mirror" . In- 
deed, the ideal mirror must support a surface current 
J„i oc eNcjac^ /to growing with the 7— factor. A realistic 
plasma surface does not provide such a current and the 
boundary condition must be changed. 

Let us consider the tangential vector potential compo- 
nent of a laser pulse normally incident onto an overdense 
plasma slab: 



1 d^A{t,x) d^Ait,x) 
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where A{t, x — —00) = and j is the tangential compo- 
nent of plasma current density. Eq. yields 



-t-00 

A{t,x)^2n J 3 {t,x,t',x') dt'dx'. 



(2) 



Here 3{t,x,t',x') = j (i', a;') (e_ - 6+). We have 
defined 6- = e {t - t' - \x - x'\/c) and 6+ = 
Q {t — t' + {x — x')/c), ^{t) is the Heaviside step- 
function. Due to this choice of J the vector potential 
A(t, x) satisfies both Eq. and the boundary condi- 
tion at X — —00 since J{t,x = —oo,t',x') = 0. The 
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FIG. 1: A sketch of possible motions of the plasma surface, 
(a) Single "hump" per half-laser period: two saddle points 
are responsible for the generation of nth harmonic, (b) More 
complex plasma surface motion, the saddle points appear in 
pairs, two per "hump". Interference between the different 
saddle points leads to modulations of the spectrum 0. 



tangential electric field is Ej = —{l/c)dtA{t,x). At the 
electron fluid surface X{t) we have 



1 



j{t + a^/c,X{t) + OdC (3) 



where £^ = x' — X{t). If the characteristic time r of 
the skin layer evolution is long in the sense ct ^ S, 
where d is the plasma skin length, then we can Taylor- 
expand: j{t ± £^/c,x' = X{t) + S^) « i{t,x') ± e, where 
e = {^/c)dtj{t,x'). We substitute this expression into 
(0. The zero order terms cancel one another and we 
get Et{t,X{t)) oc Jp{S/cT) El, where Jp oc cEi is the 
maximum plasma surface current. Thus, as long as the 
skin-layer is thin and the plasma surface current is lim- 
ited, we can use the Leontovich boundary condition 



Et(i,X(t)) = 0. 



(4) 



The same boundary condition was postulated ad hoc by 
Plaja et al^ to interpret their PIC simulation results. In 
the present work, we substantiate the boundary condition 
Q physically. Using it we are able to derive very general 
properties of the reflected radiation. 

According to Eq. Q, the reflected wave electric field 
at the plasma surface is F,r{t' , X{t')) = -E,{t' ,X{t')), 
where F,^{t' , X {t')) = -{l/c)dt'A,{t' , X{t')) is the in- 
cident laser field, t' is the reflection time. The one- 
dimensional (Id) wave equation simply translates a signal 
in vacuum. Thus the reflected wave fleld at the observer 
position X and time t is Er(t, x) = —Ei{t', X{t')). Setting 
a; = at the observer position we find that the Fourier 
spectrum of the electric field Er(t, x — Q) coincides with 
the spectrum of F{t) = (Aquj/c) cos{2ujt' — Lot), where 
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FIG. 2: Spectra of the reflected radiation for the laser am- 
plitudes ao = 5, 10, 20. The broken line marks the universal 
scaling I o^ . 



is the retardation relation. The fine structure of the 
spectrum of F[t) depends on a particular surface motion 
X{t), which is defined by complex laser-plasma inter- 
actions at the plasma surface. The previous theoretical 
works on high order harmonic generation from plasma 
surfaces 0,13 tried to approximate the function X{t). 

It appears, however, that a universal spectrum scal- 
ing can be obtained without an exact knowledge of the 
function X{t). Being interested in universal results we 
discard calculating X{t) and this makes our approach 
very different from the previous ones. We only sup- 
pose for a moment that the boundary motion is pe- 
riodic X{t + tt/lo) — X{t). Later we consider non- 
monochromatic laser pulses and get rid even of this re- 
striction. First, we mention that the Fourier spectrum of 
F{t) can be represented in the form: 



Fn^AoUJ Cn+l + Cn-l - i{Sn+l - Sn-l) /(2c), (6) 

where C„ and S'„ are the n— th harmonics of C{t) = 
cos(2(jj/:') and S{t) = sin(2wi'), and t' is the retarded 
time We examine only the spectrum of C{t), because 
that of S{t) can be worked out analogously. 

It is easy to see that the function C{t) has a period 
tt/w. Thus, its spectrum contains only even laser har- 
monics. Introducing the new variable t = 2t we obtain 



[exp (i$i(r)) + exp (i$2(r))] dr, (7) 



t' - X{t')/c^t. 



(5) 



where $i(r) = {l-m)LJT+2^r{T), $2(7") = -(l+m)cjT- 
2<I>r(T), and ^^(t) = {uj/c)X (arccos C (r/2) /(2a)))is the 
retarded phase. The definition of ^rir) is recurrent, be- 
cause C{t) itself is defined through X. 

To calculate the spectrum {Tj) we use the saddle point 
technique. The saddle points r„ are obtained from the 
equations (i<I>i(r„)/dT = and d^2{Tn)/dT = 0. The 
first equation reduces to d^r{Tn) / dr = n/2 — 1. Taking 
the first derivative of the recurrent definition of $r we 
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re-write the saddle-point equation as 

dX{Tr,)/dT = c{l-2/n) or n « 472(T„), (8) 

where T(t) = (l/2w) arccosC(T/2), r„ = r(r„) and 
7 = 1/^1- /32^ c/3 = dX{T)/dT. Eq. © has a clear 
physical meaning. The reflected radiation frequency is 
multiplied by the factor 47^ because of the relativistic 
Doppler effect, where 7 is the relativistic factor of the 
plasma surface. If the plasma surface oscillates non- 
relativistically, so that 7 « 1, then Eq. © has no real 
solutions for n > 1, and the spectrum of C{t) exponen- 
tially decays. When 7 3> 1, there is a real solution for 
any n < ric = ^I'f-nay.- ^ similar consideration shows that 
the second saddle-point equation has no real solutions. 
The spectrum is thus 

^n^Y. v/27r/M?*r(r„)| exp (z$i(t„) ± i7r/4) . (9) 

To estimate ^ we use the fact that the highest har- 
monics are generated around the time Tmax, when the 
plasma surface moves toward the laser with the highest 
relativistic factor 7max- In its vicinity, one can approxi- 
mate 7 « 7max[l -g^{T- Tmax) V2] as shown in Fig.Uli. 
A straightforward algebra leads to 



\dt^r{Tn)\ = Dgujn'^y^n{nc - n)/nc + 0{n^), (10) 

where D = sn\{2LoTmax) I'^lrnax- The estimate (|10|) leads 
to the spectrum intensity scaling 

Cl cx for n < n^. (11) 

When one considers the physical mechanism of high- 
harmonic generation as presented in Fig. ^ it becomes 
evident that the harmonics are emitted in the form of 
subattosecond pulses. Indeed, all harmonics above a 
number n < ric are generated at times between Tn 
and Tn'^K From Eqs. ©, (|HJl the subattosecond pulse 
duration At = t(ri^-') — t{Tn~^) as it is seen by the ob- 
server is At = (A.lZg'n)^/ Ticln for 1 <C n <C n^- This 
estimation tells us that the reflected pulse can be made 
very short by applying a filter selecting harmonics with 
high numbers larger than n. 

The mechanism presented in Fig. ^ has another very 
interesting consequence. Each harmonic is generated 
due to the saddle points corresponding to the proper 7- 
factor. These saddle points come into © with different 
phase multipliers. Fig.^ represents the case of a more 
complicated plasma surface motion, when the 7— fac- 
tor has several maxima, "humps" , per half-laser period. 
Yet, one can see that the saddle points are grouped, i.e. 
there is a couple of the saddle points on every " hump" . 




FIG. 3: Electron distribution function The helix represents 
the electron surface motion in the laser field. The reddish 
downward spikes stay for the surface relativistic motion to- 
wards the laser. These spikes are responsible for the zeptosec- 
ond pulse generation. 

The phase shift between the contributions of the sad- 
dle points t\^^ belonging to the the same ith "hump" is 

not very large. For example, for the saddle points T^^' 
this phase shift equals ^(n) = —(^LojZg^^Jn^n. The 
frequency modulation period due to their interference is 

(X 27rtj/(d$(ri)/dn), i.e. VL oc (67rg/5^7i;)n3/2. This 
inteference really brings modulation to the spectrum only 
if ^ QTTg/{bLjj). On the other side, the phase shift 
between the contributions from different "humps" can 
be much larger. This means that a non-trivial motion 
of the critical surface producing more than one 7-factor 
"hump" per oscillation period is the cause of the spec- 
trum modulation. This conclusion supports the explana- 
tion of the modulation proposed by Watts et aZQ and 
agrees with the experimental observations by Teubner 
et id () and Watts et aZQ- Finally, we notice that the 
larger number of the saddle points does not change the 
averaged value for c?($(i„) and, consequently, does not 
affect the overall spectrum scaling. 

A careful analysis reveals that the intensity spectrum 
oc valid for a monochromatic incident wave 

only. Indeed, if the laser pulse is short and has a spec- 
tral bandwidth Aw, then the spectral scaling fails at 
frequencies > uj'^/Auj. At these frequencies, the re- 
flected radiation is broadband and contains no distin- 
guishable harmonics. To examine the finite bandwidth 
influence we represent the incident wave as Ei(i,x) = 
Re [ujAq eKp{—i(f)(t + x/c))], where dt4>{t + x/c) w uj and 
Aoj <C A procedure almost identical to that lead- 
ing to Eqs. ITJ) and © reduces the Fourier spectrum 
calculation of E,. to that of C{t) = cxp {—i(f> {T{t))), 
where T{t) = -t + 2X{t')/c, 2t' ^ t + F (i InC(t)) and 
F{<j){t)) = t for any t. 

To apply the steepest descent method to evaluate the 
il-Fourier amplitude Cq we introduce a new function 
$ = Re(0) such that dt$(-t + 2X{t')/c) = n. We have 
neglected Im((/)) because the shift of the saddle point due 
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FIG. 4; Zeptosecond pulse train: a) temporal structure of 
the reflected radiation; b) zeptosecond pulse train seen after 
spectral filtering; c) one of the zeptosecond pulses zoomed, its 
FWHM duration is about 300 zs. 



to lm{(j>) is only 0{{uj/fl)'^). Instead of the monochro- 
matic expression © we obtain now (x 1/ i/ldj $|, 
where 



B + O 



guj / {2ujnc — ^) 



In 



(12) 



for VI < 2u!nc. Here, the 0{. . . )-term on the right 
hand part of Ea. (|12(l is inherited from the monochro- 
matic approximation, see Eq. IjlOjl . The new term B — 
\2ioC^axd/dT [{d^{T,r,ax)/dT)/{d(^{T^ax)/dT)]\ appears 
due to the finite spectral width. Here Tmax, and g are 
the same as in Eq. ((Tn|l . and Cmax = C{Tmax)- 

One readily finds that B (x l/r^, where r oc l/Aw is 
the incident pulse duration. Thus, we can neglect the 
0{. . . )-term for Q/lu > rid — (gt^) / (jmax^^) what leads 
to 



In oc 



Cn 



(13) 



We emphasize that the criterion Q > rid defining the 
transition from l/w^^^-spectrum to l/cj'^-spectrum de- 
pends on the function B that will be analised in detail 
elsewhere. Altogether, the spectral intensity /„ of the 
reflected radiation scales as In oc l/n^/"^ for 1 < n < rid, 
then as Jo oc 1 /n^ for rid < n < ^'^f^^^^^ and finally decays 
exponentially for n > Aj^ 



To check our analytical results, we have done a number 
of Id PIC simulations. A laser pulse with the Gaussian 
profile a = oq exp {—t^/rf) was incident onto a plasma 
layer with a step density profile. Fig. |21 shows spectra of 
the reflected radiation for laser amplitudes ao — 5, 10, 20, 
duration ljtl — 47r and the plasma density N,, — 307Vc, 
which roughly corresponds to the solid hydrogen or liquid 
helium. The log-log scale of Fig. [3 reveals the power-law 
scaling of the spectral intensity In oc 1/n^. The crit- 
ical harmonic number ric, where the power-law scaling 
changes into the exponential decay increases for higher 
laser amplitudes ._ Also, the spectral intensity modula- 
tions are seen 



/ max 



Let us take a close look at the particular case oq — 20 
(the red line in Fig. I^J. In this case, the power-law spec- 
trum extends above the harmonic number 2000, and zep- 
tosecond pulses can be generated. As one sees from the 
electron distribution function fit, x,px), Fig.|2| the max- 
imum surface 7— factor 7inax ~ 25 is achieved at the time 
t « 6. The temporal profile of the reflected radiation 
is shown in Fig. 0] When no spectral filter is applied. 
Fig. a train of attosecond pulses is observed . How- 
ever, when we apply a spectral filter selecting harmonics 
above n ~ 300, a train of much shorter pulses is recov- 
ered. Fig. Fig. ^ zooms to one of these pulses. Its 
full width at half maximum is about 300 zs. At the same 
time its intensity normalized to the laser frequency is 
huge [eEzs/mcLoY ~ 14 that would correspond to the 
intensity /^^ « 2 x lO^^ W/cm^. 

The presented theory and simulations are Id, i.e., in a 
planar geometry. This assumes that the laser focal spot 
radius r ^ A. Recent 2d and 3d PIC simulations in 
the regime j3] have shown that the multi-dimensional 
effects may help, particularly, to isolate a single (sub- 
) attosecond pulse. 

In conclusion, the reflection from ultrarelativistically 
moving plasma boundaries forms two different universal 
spectra because of the relativistic Doppler effect: 1/n^/^ 
for an incident quasimonochromatic wave and l/rv^ for 
a broadband pulse. The slow power law decay of the 
high harmonic spectrum causes high intensity zeptosec- 
ond pulses. The observed pulse duration of 300 zs corre- 
sponds to the spatial pulse extent comparable to the first 
Bohr orbit. As a result, the special relativity theory in- 
trodicing the relativistic Doppler effect may open a new 
page of applied ultrafast laser spectroscopy. 
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